Anharmonic properties of double giant dipole resonance 
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A systematic microscopic study of the anharmonic properties of the double giant dipole resonance 
(DGDR) has been carried out, for the first time, for nuclei with mass number A spanning the whole 
mass table. It is concluded that the corrections of the energy centroid of the = + and 2 + 
components of the DGDR from its harmonic limit are negative, have a value of the order of few 
hundred keV and follow an A -1 dependence. 
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Many-fermion systems, from metals in bulk to atomic 
nuclei, display collective degrees of freedom known as 
plasmons in the language of solid state physics [Q and as 
giant resonances within the framework of nuclear physics 
P). A very successful description of these modes is pro- 
vided by the time-dependent mean-field theory in 
the variety of versions known as quasi-boson approxima- 
tion, random phase approximation, linear response the- 
ory , time-dependent Hartree-Fock theory , time- 
dependent local density approximation [fjj, etc. At the 
basis of all these, to a large extent, equivalent theoretical 
descriptions of the sloshing back and forth of electrons 
against ions, of protons against neutrons, etc, one finds 
the small amplitude approximation which identifies these 
modes with the one-phonon state of an harmonic oscilla- 
tor. While it is true that a couple of fermions (particle- 
hole excitation) has integer spin, its behavior cannot be 
identified for all relative energies and momenta with that 
of a real boson, the range of validity of this identifica- 
tion depending on the correlation energy of the pair. 
On the other hand, close to the ground state, fermion 
particle-hole excitations do behave as (quasi)bosons. In 
fact, the terms which, in the equations of motion, are 
related to the non-bosonic contributions of the commu- 
tation relations of pairs of fermions have random phases 
leading to cancellations which reduce conspicuously the 
contribution of the corresponding terms, eventually jus- 
tifying the harmonic approximation In any case, 
all degrees of freedom of a many-fermion system are ex- 
hausted by the degrees of freedom of the particles. Conse- 
quently, although collective vibrations display small over- 
laps with each of the (particle-hole) components of the 
wave function describing the mode, a certain amount of 
over-counting is unavoidable. 

With the advent of high fluency lasers and of high lumi- 
nosity heavy ion beams at relativistic energies, it is now 
possible to study multi-plasmon states in bulk matter 
and in clusters [|f0| , as well as states of multiple excited 
giant dipole resonances in atomic nuclei pd] , and thus 
test the limits of validity of the harmonic paradigm in 



many-fermion systems fl2|| . In particular, the discovery 
of the double giant dipole resonance in nuclei (DGDR) 
]iT| , ^8| -[2l]] and the observation of small deviations from 
the harmonic picture concerning the excitation energy 
and the spreading width, combined with the large (up 
to a factor 2-3 enhancement) deviations of the associ- 
ated Coulomb excitation cross sections measured in rel- 
ativistic heavy ion collisions (f^| , call for a better under- 
standing of the role anharmonicities play in the spectrum 
of the DGDR. In fact, anharmonicities influence electro- 
magnetic DGDR cross sections in several ways: a) the 
energy shifts of the DGDR states from the harmonic val- 
ues can affect in an important way the electromagnetic 
cross section, in keeping with the exponential dependence 
of these quantities with the Q- value of the process [0 , 
b) anharmonicities lead to changes in the El transition 
matrix elements to preserve the energy weighted sum rule 
(EWSR) |^3| which eventually reinforce these effects, c) 
anharmonicities which are a consequence of the mixing 
of states with different number of phonons, give rise to 
many paths, other than the (harmonic) two-step one, to 
excite the DGDR in electromagnetic processes. While all 
these questions inspired much theoretical work, [p4| [34[ , 
no clear picture has emerged of the DGDR anharmonic- 
ity question, let alone an explanation of the "Coulomb 
excitation anomaly" . In particular, no consensus exists 
concerning the mass-number dependence of the energy 
shifts from the harmonic values. 

In this Letter we present the results of the first, sys- 
tematic calculation of the spectrum of the DGDR, carried 
out in a complete one- and two-phonon basis (the effect of 
the 3 phonon states on the anharmonicity being arguably 
small p8| ) for nuclei with mass number A spanning the 
whole mass table. It will be concluded that the energy 
shift (lowering) of the energy centroid of the J n = + and 
2 + components of the DGDR from the harmonic limit is 
rather modest (few hundreds of keV) and display a clear 
A^ 1 dependence. The solution of the Coulomb excitation 
anomaly is thus likely to be found elsewhere |35| . 

The Hamiltonian used in describing the system con- 
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tains, aside from a mean field term which determines 
the single-particle motion of protons and neutrons, a 
monopole pairing interaction and a separable multipole- 
multipole force with strengths adjusted so as to repro- 
duce the odd-even mass differences and the spectrum of 
low-lying vibrations and of giant resonances respectively 
P,[7| JT5|Jl7| | . In particular, the strength of the isovector 
dipole-dipole term was fixed by fitting the observed en- 
ergy centroid of the GDR in each nucleus or, lacking this 
information, the value emerging from the energy system- 
atics (80 -A- 1 / 3 MeV). 

The basis of one-phonon states was obtained by diag- 
onalizing the Hamiltonian in the quasiparticle random 
phase approximation. The basis of two-phonon states 
was constructed by coupling two one-phonon states to 
total angular momentum and parity J n = + and 2 + , 
in keeping with the quantum numbers of the DGDR 
states. The two-phonon basis thus include, aside from 
the states [1~ x 1jv]o+(2+)i where the subindcx i is used 
to distinguish between the different one-phonon dipole 
states arising from the shell structure of the system, also 
two-phonon states made up of + , 1 — , 2 + , 3~, and 4 + 
phonons. All one-phonon states, up to 40-50 MeV of ex- 
citation energy and contributing with more than 1.0% to 
the EWSR (0.2% in the case of dipole modes) have been 
included in the calculations. This choice leads, for the 
J w = 2 + component of the DGDR in heavy nuclei, to a 
two-phonon basis containing of the order of 10 3 states. 

The Hamiltonian written in terms of quasiparticles and 
phonons jlTj] is diagonal in the space of one- and two- 
phonon states separately, but contain terms coupling one- 
to two-phonon states. Diagonalizing this Hamiltonian, 
we obtain the total wave functions \&j and the corre- 
sponding eigenvalues from which the results displayed in 
Table Q have been obtained. In the second column of this 
table, the percentage of the Thomas-Reiche-Kuhn EWSR 
exhausted by the selected one-phonon dipole states is dis- 
played, while the third column contains the percentage 
of the EWSR for the DGDR (calculated in Ref. ||) 
exhausted relatively to the sum of the + and 2 + compo- 
nents. The small differences observed between the per- 
centage of the EWSR exhausted by the DGDR and the 
GDR is mainly due to the fact that the ground state 
is considered in the calculations as the one-phonon vac- 
uum, the ground state correlations arising from the in- 
teraction between multi-phonon configurations not being 
taken into account. In the columns four and five the en- 
ergy shifts AE c (J n ) of the centroids of the J n — + and 
2 + members of the DGDR with respect to the harmonic 
predictions are reported. In Fig. |l| we show the quantity 



+ and 2 + of the nucleus 136 Xe. The same calculations 
have been repeated in a basis containing a single two- 
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carrying the largest fraction of the EWSR, and all one- 
phonon states so as to reproduces as far as possible the 
harmonic scenario within the framework of the present 
microscopic calculation. We shall discuss these results 
before discussing those of the full calculation. 

The diagonalization in the reduced space leads to a 
breaking of the [1~ x l~]o+(2+), and thus to a set of 
states with J 71 " = + and 2 + , one of which carries about 
95% of the two-phonon configuration oscillator strength. 
The energy shift of this state from the energy of the 
(non-interacting) two-phonon configuration is shown in 
Table |j] in the columns labeled "Sum" . There are two 
mechanisms contributing to this shift: the first one is 
associated with the Pauli principle corrections. Exclud- 
ing four-quasiparticle configurations which violate Pauli 
principle reduces somehow the collectivity of two-phonon 
configurations. One thus expects a downward shift for 
isovector phonons like e.g. the DGDR (cf. column I of 
Table displaying the results obtained including only 
Pauli principle like-processes). This shift is found to 
scale with A -1 as expected from general arguments |3(| 
and simple models Ea,p3 34 39] . The second mechanism 
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tor the different (two-phonon) states is, eigenstates of the 
total Hamiltonian with angular momentum and parity 



arise from the interaction of the [1 
tion with all one-phonon states. The energy shifts arising 
from this interaction are given in columns II of Table ||. 
A strong cancellation with the first contribution is found, 
although not as complete as that reported in ref. p3[ , 
where estimates of the two contributions to the total en- 
ergy shift under discussion have been carried out within 
a schematic model. This (second) contribution arising 
from the interaction of two-phonon configurations with 
one-phonon states is found, in the present simplified cal- 
culations, not to have any simple dependence with A, as 
it arises from the coupling of the single two-phonon con- 
figuration chosen, to relatively few one-phonon configu- 
rations lying close in energy and displaying a moderate 
value of the coupling matrix elements. 

Carrying the diagonalization in the full two- and one- 
phonon space, the contribution to the DGDR energy cen- 
troid associated with the second mechanism vanishes be- 
cause the trace of the matrices is conserved independently 
on the values of non-diagonal matrix elements. Conse- 
quently, the centroid of each of the DGDR configurations 
(diagonal elements) remain at the non-interactive value, 
eventually modified by the Pauli principle for two-phonon 
corrections. The corresponding energy shifts AE C (J V ) 
for the double magic nuclei 40 Ca and 208 Pb obtained 
in the present calculation (cf. Table | columns four 
and five) are of the same order of magnitude as those 
obtained in microscopic calculations in Refs. |24|,|26||52]] . 
On the other hand, the mixing between one- and two- 
phonon states obtained for the other nuclei (cf. Ta- 
ble [n]) are larger than for 40 C a and 208 Pb. This is in 
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keeping with the fact that doubly-magic nuclei are much 
more rigid than the semi-magic ones. The most collec- 
tive [17" x l, r ]o+(2+) configuration in general prefers to 
mix with either [LEORxHEOR] or other [17" x l i 7]o+(2+) 
configurations where L(H)EOR is the low (high) energy 
octupole resonance. 

From the systematic calculations one can extract the 
A-dependence of the energy shifts of the centroids from 
the harmonic limits. For this purpose, the value of the 
energy shifts of the DGDR are shown in Fig. || as a func- 
tion of A. The continuous curve represents an A^ 1 fitting 
to the data, while the dashed line shows the A~ 5 ^ 3 de- 
pendence obtained in the variational time-dependent ap- 
proach of Ref. |3lJ] . The results of our calculations follow 
quite accurately the A -1 behavior, even if both doubly- 
and semi-magic nuclei have been included in the system- 
atics. Weighting equally the + and 2 + components of 
the DGDR we obtain from a \ 2 analysis of the results 
displayed in Fig. ||, AE = b ■ A~ a with a = 1.08 ± 0.06 
and b = -37 ± 8 MeV. 

Although our conclusion on the A -1 dependence is 
based on calculations within a specific model, general ar- 
guments |3(| and estimates [ ^5| , ^3||3l|j3^| support it. Dif- 
ferent ^4-dependence of the energy shift are discussed in 
fiijl , in terms of the number of active nucleons. The 
present results indicate that in the case of the GDR this 
number is indeed of the order of A 2 / 3 , as for the f2 factor 
in ||. 

We conclude that the deviation of the energy centroid 
of the double giant dipole resonance from the harmonic 
limit displays a behavior with mass number A typical of 
that associated with the global properties characterizing 
the system, like e.g. the energy centroid of the giant 
dipole resonance. 

Discussions with G. F. Bertsch, K. Hagino, I. Hama- 
moto, and B. Mottelson are gratefully acknowledged. We 
thank J. Bryssinck for help. V. Yu. P. acknowledges a 
support from INFN and NATO. 
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FIG. 1. Energy distributions of the B(E1 x El) values as- 
sociated with the excitation of the + and 2 + components of 
the DGDR in 136 Xe, in comparison with the same quantity 
for the 2 + component in the harmonic limit. Scales are chosen 
proportionally to (2 J + 1). 
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FIG. 2. Shift of the DGDR energy centroid (0 + - stars and 
2 + triangles) from the harmonic limit. The continuous and 
dashed curves represent fits, assuming an A -1 and an A~ 5 ^ 3 
dependence respectively, of the results of the microscopic cal- 
culations. 

TABLE I. Percentage of the EWSR exhausted by the GDR 
and DGDR of the atomic nuclei indicated in the first col- 
umn. In column 4 and 5 is displayed the anharmonicity shift 
AE C (J W ) of the energy centroid of the J w = + and 2 + com- 
ponents of the DGDR from its harmonic limit. 
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TABLE II. Energy shift of the two states [lr x ir] J7r 
(J* — + and 2 + ) with respect to the harmonic value 
27iaj(l~). The label io indicates the component of the GDR 
carrying the largest fraction of the EWSR. The calculations 
have been carried out in a basis which includes only the 
two-phonon configuration [lr x l^jj" and a complete set of 
+ (2 + ) one-phonon states. The contributions to the energy 
shift arising from Pauli principle corrections and due to the 
interaction of the two-phonon configuration with one-phonon 
configurations are shown separately in I and II, respectively. 
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